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^t ■ Abstract 

G\ 

The connection is established between two different action principles for 
perfect fluids in the context of general relativity. For one of these actions, 
S, the fluid four-velocity is expressed as a sum of products of scalar fields 
and their gradients (the velocity-potential representation). For the other 
action, S, the fluid four-velocity is proportional to the totally antisymmetric 
product of gradients of the fluid Lagrangian coordinates. The relationship 
between S and S is established by expressing S in Hamiltonian form and 
identifying certain canonical coordinates as ignorable. Elimination of these 
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f*"*- ■ coordinates and their conjugates yields the action S. The key step in the 

^-r ■ analysis is a point canonical transformation in which all tensor fields on space 

Q*>\ . are expressed in terms of the Lagrangian coordinate system supplied by the 

q | fluid. The canonical transformation is of interest in its own right. It can be 

applied to any physical system that includes a material medium described 

by Lagrangian coordinates. The result is a Hamiltonian description of the 

system in which the momentum constraint is trivial. 
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I. INTRODUCTION 

Many different variational principles have been devised that describe perfect fluids in 
the context of general relativity. Some make extensive use of the Lagrangian coordinate 
system, the system of spacetime coordinates defined by the fluid flow lines and the proper 
time along those flow lines Jl| . The use of special coordinates violates the spirit of general 
relativity More importantly, the Lagrangian time coordinate does not, in general, label 
spacelike hypersurfaces and is not suitable for the development of a Hamiltonian version 
of the action. Some perfect fluid actions rely on constrained variations in which the field 
variables are not freely varied |2|. In practice it is unclear whether or not such an action can 
be expressed in Hamiltonian form. Moreover, the usual relationship between the variational 
principle and its associated boundary value problem || is absent because typically the 
boundary conditions that follow from the action must be imposed on the parameters of the 
constrained field variations rather than on the field variables themselves. 

Some perfect fluid actions employ a set of spacetime scalar fields as the basic field vari- 
ables. The fluid four-velocity is then constructed from combinations of these fields and their 
gradients. This is the velocity-potential or Clebsch potential [ffl formalism — for general rel- 
ativists perfect fluids this approach was first discussed by Tarn and O'Hanlon [|| and by 
Schutz 0. (See also Ref. 0.) In these papers the theory of Pfaff forms || was used to 
minimize the number of velocity potentials. Thus, the fluid four-velocity U a was expressed 
in terms of four scalar fields, <p, #, W, and Z, as 

fiU a = -(p^-s^ + WZ ta , (1.1) 

where /i is the fluid chemical potential and s is the entropy per particle. The velocity 
potential formalism was further developed in Ref. where it was argued that locally the 
fluid velocity can be expressed as in Eq. (1.1), but globally the representation (1.1) must 
allow for potentials that are not single-valued functions on spacelike hypersurfaces. On the 
other hand, if the velocity potentials (p, •&, W, and Z are required to be single-valued then 
the representation (1.1) places a restriction on the velocity U a . In particular, it can be 
shown that the helicityQ for an isentropic fluid (s = const) necessarily vanishes when the 
velocity is given by expression (1.1) with single-valued potentials. 

The shortcomings of the velocity-potential representation (1.1) are avoided by enlarging 
the number of potentials and writing |J 

pU a = -<p >a - s$, a + W k Z\ a , (1.2) 

where the index k takes the values 1, 2, 3. Thus there are eight scalar potentials </?, $, 
Wk, and Z h . In this case the fields Z k can be interpreted as Lagrangian coordinates for the 
fluid. That is, the values of the three fields Z k {y) serve as a set of labels that specify which 
fluid flow line passes through the spacetime point y. Perfect fluid actions that are based on 



1 Helicity is defined as the integral over a spacelike hypersurface of —V A dV, where V := {J.U is 
the Taub current [|i| expressed as a one-form. 



the velocity-potential representation (1.2), with Z k interpreted as Lagrangian coordinates, 
will be denoted by S. There are several such actions corresponding to different functional 
expressions for the fluid equation of state ||. These actions have several features in common 
with yet another class of perfect fluid actions first discussed by Carter jll] and by Kijowski 



et al. |I2| . (See also Ref. pif .) These actions will be denoted by S. They are distinguished 
by the fact that they are functionals only of the Lagrangian coordinates Z k . In this case, the 
four-velocity is constructed from the antisymmetric product of gradients of the Lagrangian 
coordinates, 

jja _ e ^P Z \pZ\ a Z 3 tP , (1.3) 

where e af3crp is obtained by raising indices on the spacetime volume form e a p ap . 

The results obtained in Ref. [(J suggest that the actions 5* and 5* are closely related. The 
first objective of this paper is to establish this relationship. 

In brief, the connection between the actions S and S is described as follows. Both of 
the actions S and S can be put into standard Hamiltonian form in which the Hamiltonian 
is a linear combination of contributions to the Hamiltonian and momentum constraints. 
For S, the variables s and Wk are proportional to the momenta conjugate to $ and Z k , 
respectively, so the phase space variables are ip, $, Z k , and their canonical conjugates. By a 
point canonical transformation, the Hamiltonian action S can be expressed in terms of a new 
set of canonical variables. The new canonical coordinates <p and i9, that in an appropriate 
sense replace (p and ■&, are ignorable (cyclic) and can be eliminated along with their canonical 
conjugates by Routh's procedure |J. The remaining variables, after application of the inverse 
canonical transformation, are the Lagrangian coordinates Z k and their conjugates. The 
resulting action is the Hamiltonian form of S. 

The key step in the analysis is the point canonical transformation. This transformation 
is described in terms of a mapping Z : S — > S from the space manifold £ to the abstract 



"fluid space" manifold S whose points ( E S are the individual fluid flow lines [|12j,p||. The 
canonical variables Z k (x) collectively express the mapping Z in terms of local coordinates 
x a on £ and ( k on S: 

Z:Y,^S by x a i-> ( k = Z k {x) . (1.4) 

Under the canonical transformation, the canonical coordinate Z k {x) is replaced by its inverse 
X a ((), which is the local coordinate expression of the inverse mapping X := Z~ x . Whereas 
(-k _ ^(^x) specifies which flow line passes through the space point x a , x a = X a (() specifies 
the spatial location of the flow line ( k . 

Under the canonical transformation all other tensor fields are mapped from space E to 
the fluid space S by Z and its inverse X. Accordingly, the original set of canonical variables 
are referred to as the S-variables, while the new set of canonical variables are referred to 
as the iS-variables. In particular, the S-variables f(x) and i9(x) are mapped into the <S- 
variables <p(() '■— (p(X(()) and fl(() := $(X(()). Their canonical conjugates IL V (C) and 
n,?(£) are the fluid particle number per unit coordinate cell d 3 ( and the fluid entropy per 
unit coordinate cell d 3 (. Since the particle number and entropy do not change within a 
flow tube (defined by a neighborhood of a point ( in the fluid space), the momenta II ¥ ,(^) 



and II#(£) are conserved. In the variational principle the conservation of 11^ and II# is a 
consequence of the ignorable property of (p and $. 

The second objective of this paper is to present the details of the canonical transformation 
just described. This transformation can be applied to any system that includes a material 
medium described by Lagrangian coordinates Z k . The result is a canonical formulation of 
the system in which the variable X a (() and its conjugate P a {C) are the only ones that depend 
on the space manifold S for their definitions. Consequently the momentum constraint, which 
is the canonical generator of spatial diffeomorphisms, depends only on the variables X a and 
P a . Thus the momentum constraint expressed in the canonical S- variables is relatively 
trivial. This result, expressed in a slightly different form, was used in Ref. |14| to simplify 
the formal Dirac quantization of gravity coupled to a pressureless perfect fluid (dust). In 
particular, the operator condition on the quantum state functional that follows from the 
classical momentum constraint is trivially solvable. The canonical transformation between 
E-variables and iS-variables is also utilized in Ref. 111511, where the action for an elastic 



medium of clocks [|1(J is analyzed. 

Two perfect fluid actions of the velocity-potential type are presented in Section 2a. 
For one of these actions, the equation of state is specified by giving p(n, s), the fluid energy 
density p as a function of particle number density n and entropy per particle s. For the 
other action, the equation of state is specified by giving p(p,s), the fluid pressure p as a 
function of chemical potential p and entropy per particle s. Symmetries of these actions that 
preserve the form of the velocity-potential expression (1.2) are presented. The Hamiltonian 
form of the action with p(n, s) is given in Ref. f9j, while the Hamiltonian form of the action 
with p(p, s) is derived in Appendix A. These results are merged in Section 2b. The result is 
an action S whose associated Hamiltonian is implicitly defined through the solution to an 
auxiliary equation that relates the Eulerian and Lagrangian particle number densities. The 
distinction between the choice of p(n, s) or p(p, s) as the equation of state arises through the 
choice of n or p as the independent variable in the auxiliary equation. Conserved Noether 
charges associated with the symmetries of the action are displayed. The discussion in Section 
2c addresses the results obtained from an equation of state in which n and temperature T 
(or p and nT) are the independent thermodynamical variables. 

The canonical transformation from the S-variables to the 5-variables is presented in 
Section 3a, and applied to the perfect fluid and gravitational actions in Section 3b. A key 
mathematical derivation is presented in Appendix B. The momentum constraint and its role 
as the generator of space S diffeomorphisms is discussed in Section 3c. The conserved charge 
that generates fluid space S diffeomorphisms is discussed in Section 3d. 

In Section 4a the canonical coordinates 93 and 1? are identified as ignorable variables 
and the pairs (f, H^ and 1?, II# are eliminated by a straightforward application of Routh's 
procedure. The resulting action is denoted S. The symmetries and conserved charges of 
S are discussed. In Section 4b the canonical transformation is reversed for the remaining 
variables yielding an expression for the action in terms of S-variables. The Lagrangian form 
of S is obtained in Section 4c for the equation of state p(n,s). The fluid four-velocity is 
expressed as the combination (1.3) of gradients of the Lagrangian coordinates. In Section 
4d a hybrid action is constructed by eliminating the ignorable pair tp, 11^, and retaining the 



pair i9, n# as dynamical variables. This yields the action discussed in Ref. [12 



II. PERFECT FLUID ACTION WITH VELOCITY POTENTIALS 

A. Lagrangian Actions 

The local expression of the first law of thermodynamics can be written as 

dp = fidn + nT ds , (2-1) 

where p is the energy density, n is the number density, s is the entropy per particle, p, is the 
chemical potential, and T is the temperature. Pressure is defined by p := np, — p. These 
quantities characterize the thermodynamical state of a fluid in the fluid rest frame. 

Equation (2.1) shows that the fluid equation of state can be specified by the function 
p(n, s). An action principle of the velocity-potential type that describes a relativistic perfect 
fluid with equation of state p(n, s) is 

S F iJ a ,^,$,s,W k ,Z k ]laP ] 

= Jd 4 y{- v ^p(\J\/V=l,s) + J a (<P,a + s$, a -W k Z\ a )} . (2.2) 



Here, ^ a p is the spacetime metric and the magnitude of J a is denoted by \J\ := J —J ar y a pJ t3 ■ 
The variable J a is interpreted as the densitized particle number flux vector; that is, the fluid 
four-velocity is defined by 

ja 

U a := — — (2.3) 

n V ~ 7 

where n := | J\/y/— J is the number density. The equation of motion obtained by varying 
the action with respect to J a is 

M4 = -<P,a ~ Stf, a + W k Z\ a . (2.4) 

This is the velocity-potential expression (1.2). The potential •§ that appears here is the 
"thermasy" JL7]]. The equation of motion that follows from varying s, along with the first 



law (2.1), show that the gradient of ■& along the fluid flow lines is the temperature: 

$ a JJ a = T . (2.5) 

The equations of motion obtained by varying W k imply 

Z\ a U a = , (2.6) 

so that the Lagrangian coordinates Z k are indeed constant along the fluid flow lines. From 
these results it follows that the gradient of the variable p along the flow lines is given by 

V, a U a =fi-Ts, (2.7) 

which is the chemical free energy of the fluid. 



The stress-energy-momentum tensor obtained from the action (2.2) has the form ap- 
propriate for a relativistic perfect fluid. Moreover, it is shown in Ref. || that the equations 
of motion that follow from this action correctly encode the familiar perfect fluid equations, 
namely, the conservation of particle number and entropy per particle along the fluid flow 
lines, and the Euler equation relating the fluid acceleration to the gradient of pressure. 

The first law of thermodynamics also can be written in the form 



dp = n dp — nT ds 



(2.8) 



with the energy density defined by p := np—p. This expression shows that the fluid equation 
of state can be specified by the function p(p, s) . An action principle of the velocity-potential 
type that describes a relativistic perfect fluid with equation of state p(p, s) is || 



S F [V a ,<p,'d,s,W k} Z k ;^] 



d 4 y^l{p(\V\,s) 



dp(\V\,s) 

d\V\ 



\V\-^(<p, a + s$ >a -W k Z\ a ) 



(2.9) 



The variable V a is interpreted as the Taub current ||10|| , so the fluid four-velocity is defined 
by 



U a : = 



(2.10) 



where p := \V\ is the chemical potential. The equation of motion obtained by varying V c 
yields the velocity-potential expression 



V a 



-ip >a - stf a + W k Z k yCt , 



(2.11) 



which coincides with Eq. (2.4). The velocity potentials that appear in the action (2.9) have 
the same physical interpretation as those that appear in the action (2.2). 

As shown in Ref. 0, the action (2.9) yields the correct stress-energy-momentum tensor 
and the correct perfect fluid equations of motion. A simpler version of this action is obtained 
by replacing V a with the solution (2.11) of the V a equation of motion. This yields S F = 
J d 4 y y/— 7p(|V|, s) where V a is shorthand notation for the right-hand side of Eq. (2.11). 
Both versions lead to the same action in Hamiltonian form; the derivation in Appendix A 
uses the action in Eq. (2.9) as a starting point. 

The actions (2.2) and (2.9) are invariant under certain symmetry transformations that 
preserve the form of the velocity potential expression on the right-hand sides of Eqs. (2.4), 
(2.11). These symmetries include P,|i~4| invariance with respect to a change of Lagrangian 
coordinate labels, 



Z k = S k (Z') 



k dz , k 



invariance with respect to a deformation or "tilt" of the constant-*/? surfaces, 

d&(Z) 



<P 



up + &(Z) , w' k = w k + 



dz h 



(2.12a) 



(2.12b) 
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and invariance with respect to a deformation of the constant-^ surfaces, 

& = + e{Z), W' k =W k + s?^±. (2.12c) 

The invariances (2.12a-c) give rise to conserved Noether currents and charges. For the 
invariance (2.12a), the Noether current is related to the fluid momentum in the directions 



orthogonal to cp >a + s$ ;0l fli]] . The Noether currents associated with the invariances (2.12b) 
and (2.12c) are expressions of the conserved particle number current and the conserved 
entropy current ||. 

B. Hamiltonian Action 

The Hamiltonian form of the action (2.2) is derived in Ref. 0, and the Hamiltonian form 
of the action (2.9) is derived in Appendix A. These results are summarized as follows. The 
Hamiltonian action is 

S F [<p, II V , tf, n„, Z\ P k ; g ab , N ± , N a ] 

= j dt j d 3 x(n^ + 11*0 + P k Z k - N L Hl - N a H^) , (2.13) 

where g a b is the metric on space S and N 1 - and N a are the lapse function and shift vector, 
respectively. The fluid contributions to the momentum and Hamiltonian constraints are 

H F a = n^, a + U^ >a + P k Z\ a , (2.14a) 



Hi = yyn* + HZg°*H? - ^gp . (2.14b) 

In the expression for H±, p and p are determined as functions of the canonical variables as 
follows. First, choose an equation of state, either p(n,s) or p(p,s). For the case p(n,s), 
define 

/i(n,s):= — —J— , p(n,s):=n — ^ p{n,s) ; (2.15a) 

on on 



for the case p(p, s), define 



n(,,,):=^>. (2.15b) 



Next solve the equation 



- n 



HFg ab m 
r "\ 1+ ma 2 ' 16 



for either n or p (depending on the equation of state used) as a function of 11^, H£, g a b, and 
s. Now make the identification 



s :- 



n M 



(2.17) 



for the entropy per particle s. In this way either n or /i (depending on the equation of 
state used) is determined as a function of II^, ILj, H£, and g a (,. This result along with the 
equation of state and the definitions (2.15a) or (2.15b) determine /x and p, and subsequently 
iff, as functions of 11^, ILj, H£, and g a 6. 

The momentum 11^ is the particle number on S per unit coordinate cell d 3 x. The 
auxiliary equation (2.16) relates the Eulerian number density 11^/ ^fg (the number density 
per unit proper volume) to the Lagrangian number density n. The square root factor in 



Eq. (2.16) is just the relativistic gamma factor Jl + U a g ab Ub, where U a = — if,f /(/ilL^,) are 
the spatial components of the fluid four-velocity 

The Noether charges associated with the symmetries (2.12a) are 



Q[l\ = -f d 3 xP k {x)£ k {Z{x)), i\Z):- 



dS k (Z,a) 



da 



(2.18) 



cr=0 



where H(Z, a) denotes a one-parameter subgroup of transformations (2.12a). Likewise, the 
Noether charges associated with the symmetries (2.12b) are 



Q[<f>] = f d 3 xTl v (x)cj>(Z(x)) , (j>{Z) :-- 



d<$>(Z,a) 



it, ' da 

and the charges associated with the symmetries (2.12c) are 



(2.19) 



cr=0 



Q[0] 



d 3 xU4x)0(Z(x)) , 9{Z):~- 



d@(Z,a) 



da 



(2.20) 



The charges (2.18-20) generate the infinitesimal versions of the symmetries (2.12) through 
the Poisson brackets. (The derivation of the Hamiltonian action shows that W^ = —Pk/H^.) 
In particular, the charge (2.18) is the canonical generator of DiffiS, fluid space diffeomor- 
phisms, where the fluid space vector field £ (£) is an element of the Lie algebra of DiffiS 

0- 



C. Other Equations of State 

For the sake of completeness, it is appropriate to finish this section with a comment on 
other possible expressions for the equation of state. Observe that the first law as expressed 
in Eq. (2.8) is obtained from the first law (2.1) by a Legendre transformation in n and /i. A 
Legendre transformation can be performed in T and s instead, which leads to the first law 
in the form 



d(na) = ([/, — Ts)dn — ns dT 



(2-21) 



Here a = p/n — Ts is the physical free energy. This result suggests the use of an equation 
of state of the form a(n, T). The perfect fluid action with this equation of state is obtained 
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by replacing p(n,s) in Eq. (2.2) with na(n,T) +nTs 0. The temperature T is treated as 
a new variable; the equation of motion obtained by varying T reproduces the definition of s 
from the first law, namely, 

-- 8a{n ' T) (2.22) 



dT 

The Hamiltonian form of the resulting action coincides with the Hamiltonian action (2.13) 
with na(n, T)+nTs replacing p(n, s). However, this action is not truly in Hamiltonian form, 
since the extra variable T appears. T can be eliminated by inverting Eq. (2.22) for T(n, s) 
and substituting the result back into the Hamiltonian^] Through this process the equation 
of state a(n,T) is replaced by the equation of state p(n, s), and the Hamiltonian action is 
converted back into the action (2.13) with equation of state p(n,s). In a similar manner, 
a Legendre transformation in nT and s can be applied to the first law (2.8), leading to an 
equation of state in which p + nTs is given as a function of p and nT. The corresponding 
action can be obtained from Eq. (2.9). In passing to the Hamiltonian form of this action 
the variable nT is eliminated, resulting in the action (2.13) with equation of state p(p, s). 

III. CANONICAL S VARIABLES 

A. The Canonical Transformation 

The main points of this section are best illustrated with the perfect fluid coupled to other 
fields. Thus, consider a perfect fluid coupled to the gravitational field. The action for the 
combined system is obtained by adding the gravitational action to the fluid action S F . The 
result, in Hamiltonian form, is 

S[<p, u v , <&, n„, Z k , P k , 9ab , p a \ JV- 1 , N a ] 

= I dt I d^x^cp + Uj + P k Z k +p ab g ab - N ± H ± - N a H a ) , (3.1) 

where p ab is the momentum conjugate to the metric g a b. The Hamiltonian and momentum 
constraints are 

H ± = H ( [ + H^, (3.2a) 

H a = H G a + HI , (3.2b) 

where H± and H£ are the fluid contributions (2.14) and H° and H° are the familiar gravi- 
tational field contributions. 



2 One can attempt to retain T as a canonical variable by defining its conjugate momentum Pt- 
However, this leads to a primary constraint Pt = and a secondary constraint which is Eq. (2.22) 
written in terms of the canonical variables. This pair of constraints is second class, and its elimi- 
nation is equivalent to the result obtained by solving Eq. (2.22) for T and substituting the solution 
back into the Hamiltonian. 

9 



The canonical variables that appear in the action (3.1) are the pairs Z k (x) and Pk(x), 
ip(x) and IL^x), i?(x) and H$(x), and g a b(x) &ndp ab (x). These are the canonical E-variables. 
Note that (p(x), U v (x), i9(x), H$(x), g a b( x ), &ndp ab (x) are (time t-dependent) tensor fields on 
E. The variable Z k (x) is the local coordinate chart expression of a (t-dependent) mapping 
Z : E — ► S from the space manifold E to the fluid space S. Its conjugate Pk{x) is the 
coordinate expression of a (t-dependent) mapping from E to T*S, with density weight 1 
on E. In this section an alternative set of canonical variables, the S- variables, is described. 
The iS-variables utilize the fact that on E the fluid particles naturally define a preferred set 
of spatial coordinates. 

First, introduce the t-dependent mapping X, which is the inverse of Z: 

X : S -> E , X:=Z- 1 . (3.3) 

The mapping X specifies the spatial coordinates of the fluid flow lines; that is, x a = X a (() 
is the spatial location of the fluid flow line whose Lagrangian coordinate label is ( k . The 
mappings Z and X induce mappings of tensor fields from E to S. In particular, for the 
scalars tp, d and the covariant metric tensor g ab , the corresponding tensors on S are the 
pullbacks of those fields by X: 

(3.4a) 
(3.4b) 
(3.4c) 





<p ■= x> = 


if o X 


• 




tf ;= X*$ = 


doX 






9 ■= X*g . 






se definitions become 




v(0 


= <p(x(Q) , 






0(0 


= *(X(Q) , 






9ke(0 


= x a A()x h 


AC) 9ab 


(*(0) • 



(3.5a) 
(3.5b) 
(3.5c) 

The variables flf^(C) are the components of the metric tensor 

ds 2 =g ke (()d( k d( e (3.6) 

on the fluid space S. That is, g ki (() measures the proper distance ds in space E between 
neighboring fluid particles with Lagrangian coordinate labels ( k and ( k + d( k . 

The canonical transformation from the E-variables to the iS-variables is a point trans- 
formation of the canonical coordinates defined by 

(Z k (x), <p(x), *(*), g ab (x)) ^ (X a ((), <p((), 0(C), g M (0) ■ (3-7) 

The momenta conjugate to the new coordinates X a ((), <p((), i?(C) and g k g(() are identified 
by imposing the invariance of the canonical one-form, or equivalently, the invariance of 

J d 3 x{P k (x)Z k (x) + U v (x)(p(x) + Il#(x)&(x) +p ab (x)g ab (x)} . (3.8) 

Note that x a and ( k play the role of continuous indices that label the different canonical 
coordinates and momenta. 

10 



As mappings between the space manifold £ and the fluid space manifold S, Z k [x) and 
X a (Q play a special role in the canonical transformation. Consequently a separate analysis 
will be presented for the terms P k Z k in Eq. (3.8). The remaining terms in Eq. (3.8) are 
constructed from tensor fields on E, and their analysis is presented in Appendix B. 

To begin, use the fact that Z oX is the identity on S to express Z k (x) in terms of X a (Q. 
In local coordinates, we have 



z k (x(c)) = C k 

and differentiation with respect to t yields 

~(z k Jx 



Z k (x) 



x=X(Q 



x=X(Q 



x a (() 



(3.9) 



(3.10) 



This equation can be written more compactly as Z k (X(()) = —Z k }a (X(Q)X a ((). Now 
use this result along with a change of integration variables x a = X a (() in the first term of 
expression (3.8) to obtain 



/ d 3 xP k {x)Z k {x) = f d 3 ( 

JT, JS 



dX(() 



Js 



dX(C) 



d( 



P k (x(C))z k (x(0) 
p k (x(0)z k 4x(0)x a (0 



(3.11a) 
(3.11b) 



Here, \dX (() / d(\ is the Jacobian of the change of variables.^ 

The remaining terms in Eq. (3.8) are rewritten by making use of the following result, 
which is derived in Appendix B. Consider an arbitrary ^-dependent tensor field x on ^j 
and its conjugate tensor density n. (Indices on x an d vr are suppressed, and summation of 
indices is implied in expressions such as irx-) Also let H a (x;ir, x] denote the contribution 
from x an d vr to the momentum constraint. The key result from Appendix B is 



d 3 x7i(x)x{x) 



d% t(C)x(C) " H k ((; 7T, X ] Z k , a {X{Q) X a (() 



(3.12) 



where %(C), ^(C)) an d Hk((;TT,x] are the tensors (tensor densities) on S that correspond 
to x( x )i 7r ( a; )) an d H a (x;ii, x]i respectively. 

The results (3.11) and (3.12) can be applied to expression (3.8) with the result 



3 The change of integration variables x a = X a ((), which involves the new canonical coordinate 
X a (C), is described more precisely as follows. The integrand on the left-hand side of Eq. (3.11a) 
is a three-form on S, call it ui. Its counterpart on S is the three-form u := X*uj obtained by the 
pullback mapping X*. In terms of local coordinates, uJk£m{C) = X a ^{C)X b /{C)X C ^ m {Cj ^abc{X {(,)). 
Equation (3.11a) just expresses the fundamental identity J s u = J s oj in terms of local coordinates. 
Thus, the change of variables x a = X a (Q in Eq. (3.11) arises because the integrand (the three- 
form) is being mapped from £ to S. This is the way in which the continuous index set x a is 
changed to the continuous index set C fc . 
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d 3 x{P k (x)Z k (x) + n v (x)^(x) + Tl4x)$(x)+p ab (x)g ab (x)} 

d 3 ({-H k (()Z\ a (X(())X a (0 + 1^(0^(0 + 1^(0^(0 + P H (()9ke(0} ■ (3-13) 



Here, II V ,(^), n#(£), and p M (() are the fluid space tensor densities that correspond to the 
fields n ¥ ,(x), H$(x), and p ab (x): 



n„(C) := 

P M (C) := 



0*(C) 









d( 



n v (x(0) , 
n*(x(0) , 

Z\ a {X{Q)Z\ b {X{Q)p ah {X{Q) . 



(3.14a) 
(3.14b) 
(3.14c) 



Also, -Hfc(C) is Ihe fluid space covariant vector density obtained by mapping the full mo- 
mentum constraint H„(x) from E to S: 



H fc (C) := 



5X(C) 



9C 



x%(C)tf a (x(0) . 



(3.15) 



Equation (3.13) shows that the momentum conjugate to (f(() is II (/ ,(^), which is the fluid 
particle number on S per unit coordinate cell d 3 (. Likewise the momentum conjugate to i?(C) 
is II,j(£), the fluid entropy on S per unit coordinate cell d 3 (. The momentum conjugate 
to g^eiC) * s P (0) l ne expression of the gravitational momentum in the fluid coordinate 
system ( k . Finally, the momentum conjugate to X a (Q is defined by 



P B (C) := -H k (() Z\ a (X(()) 



dX(Q 



d( 



Ha(X(C)) , 



and is proportional to the momentum constraint (3.2b) for the system. 



(3.16) 



B. The Action as a Functional of S Variables 



With the results (3.13)-(3.16), the action (3.1) becomes 
S[X a , P a , tp, n v , ■&, n#, g M , p ki , N ± , N a ] 



dt / d 3 (\ P a (()x a (0 + n v (C)^(C) + n,(C)^(C) + p^OWO 

-^(^(0)^(0 + N a (x(0)p a (0 



where the Hamiltonian constraint (3.2a) is expressed in the fluid space S as 

dX(C) 



H±(C) ■■- 



H±(X(0) , 



d( 



Hi + J^U' + Hlg^Hl - ^p . 



(3.17) 



(3.18) 
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The gravitational part H°_(Q of the Hamiltonian constraint is obtained from H G (x) by 
replacing the gravitational S-variables g a b{x), p ab {x) with the corresponding iS-variables 
9kti0i P ke (0- For the perfect fluid part of the Hamiltonian constraint (3.18), p and p are 
determined as before: First choose an equation of state p(n, s) or p(p, s) and define the 
relevant thermodynamical variables as in Eqs. (2.15). Next, solve the equation 



- n 



H F n ki H F 

vs-y+w (319) 



<p 



for n or p as a function of n^(£), H F ((), g ke ((), and s. With the identification s := II^/II^, 
then // and p are determined as functions of n v (£), n,j(£), H k ((), and g ki (()- The fluid 
space density H F (Q that appears here and explicitly in the Hamiltonian constraint (3.18) is 
the fluid contribution to the momentum constraint, mapped to S. It is expressed as a func- 
tion of the canonical 5-variables through the solution of Eq. (3.16) with the identification 
H k :=H F k +H G k . This yields 

HliO = -Pa(C)X a AC) - H G k (() , (3-20) 

where H k (() is the gravitational contribution to the momentum constraint, mapped to S. 
Thus, H G (Q depends only on the gravitational canonical variables g k e(() and 7V H ((). 

We can make a further change of variables in the action by defining fluid space tensors 
that correspond to the lapse function and shift vector: 

JV^C) := ^(1(0) , (3.21a) 

N\() := N a (X(())Z k , a (X(()) . (3.21b) 

The action then becomes 

S{X a , P a , ip, n v , 0, n*, g ki , p ke , N\ N k ] 

= f^dt f s <Pc{p a (Qx a (Q + n v (C)^(C) + n*(CW0 +p h (C)^(C) 

-N L (()H ± (() - N k (()H k (()} , (3.22) 

where the momentum constraint 

H k (():=-P a (()X a , k (()=0 (3.23) 

from Eq. (3.16) depends only on the fluid coordinates X a and their conjugates P a . 

The Hamiltonian constraint H±(x) that appears in the action (3.1) is smeared with a 
spatial scalar N^-fa). This defines a canonical generator, 

HlN 1 -] := f d 3 xN ± (x)H ± (x) = f d 3 ( ^{X^H^Q , (3.24) 

which, as shown here, can be written either in terms of the canonical ^-variables or the 
canonical iS-variables. On the other hand, the Hamiltonian constraint that appears in the 
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action (3.22) is smeared with a fluid space scalar N (£). The corresponding canonical 
generator is 

MN 1 -] := f (PC^iOHAO = f d 3 xN ± (Z{x))H ± {x) . (3.25) 

Js Jt, 

Now, the Poisson bracket of any canonical variable F with the smeared constraint HlN- 1 ] 
equals 

WHIN 1 -]} = I s £c[n\X(C)){F,H ± (C)} + {F,N ± (X(C))}H ± (C)) 

= {F,H[N ± ]} + f s <P({F,N J -(X(Q)}H ± (Q , (3.26) 

where N and N 1 - are related as in Eq. (3.21a). This result shows that the brackets 
{F, H[N ± }} and {F, iffiV- 1 ]} are equal modulo the constraints H±(x) = = H±((). Thus, 
on the constraint surface, the Poisson bracket of F with H[N } gives the change in F under 
a displacement of the hypersurface E by a proper time N (Z(x)) in the direction orthogonal 
to the hypersurface. 

C. Momentum Constraint and DiffE 

The momentum constraint H a (x) smeared with an externally prescribed vector field 
N a (x) on E can be written either in terms of the canonical E-variables or the canonical 
iS-variables: 

H[N] := J^d 3 xN a (x)H a (x) = - Jj^ N a (X(())P a (() . (3.27) 

H[N] is the generator of DiffE through the Poisson brackets. That is, for any canonical 
variable F, the Poisson bracket 

{F,H[N]} (3.28) 

gives the change in F due to an infinitesimal diffeomorphism of S generated by the vector 
field N . If F is expressed solely in terms of ^-variables, then the Poisson bracket is most 
conveniently computed using the ^-variables and the first integral in Eq. (3.27). If F is 
expressed solely in terms of iS-variables, then the Poisson bracket is most conveniently 
computed using the iS-variables and the second integral in Eq. (3.27). If F is a tensor field 
on E, then the Poisson bracket (3.28) yields the Lie derivative £ftF. 

Observe that the iS-variables <p, 11^, i9, n#, g, and p have vanishing Poisson brackets 
with H[N]. This simply reflects the fact that they are tensor fields on the fluid space S 
and are invariant under DiffE. Among the canonical iS-variables, only X a and P a have 
nontrivial transformation properties under DiffE: 

{X a ((),H[N}} = -N a (X(()), (3.29a) 

{P a ((), H[N]} = P b (()N\ a (X(()) . (3.29b) 
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The transformation (3.29a) can be verified as follows. Let M a : £ — *■ £ denote the one- 
parameter family of diffeomorphisms of E with generator 



JV: ^ 



da 



cr=0 



(3.30) 



Also let J\f* generically denote the action of J\f a on the canonical field variables. For the fields 
Z k (x), viewed as scalars on E, the action of M a is given by the pullback M*Z k = Z k oJ\f ff . 
The associated infinitesimal diffeomorphism is then 

^(KZ\x))l =o = Z\ a {x)N«{x) = {Z k (x),H[N}} . (3.31) 

Correspondingly, the action of DiffE on the mapping Z is defined by M*Z = Z oJ\f ff . Since 
X is the inverse of Z, the action of DiffE on X is given by M*X = M~ Y o X. For the 
associated infinitesimal diffeomorphism, this yields 

^(KX a (())l =0 = -N«(X(()) , (3.32) 

which agrees with Eq. (3.29a). 

The transformation (3.29b) can be checked by recalling that P a (() is a mapping from S 
to T*£ with density weight 1 on S. Thus, P a (()X a t k{() is a covariant vector density on S 
and is invariant under DiffE. Given the transformation (3.29a) for X a ((), Eq. (3.29b) is the 
required result for the vanishing of the Poisson bracket {P a (QX a t k{C), H[N]}. 

The smeared momentum constraint (3.27) appears as a term in the action functionals 
(3.1) and (3.17). The action (3.22) employs the Lagrange multiplier N (£) and defines a 
smeared constraint 

H[N] := J s d 3 (N k (()H k (() = J^d z xN k {Z{x))X\ k {Z{x))H a {x) . (3.33) 

— * — * 

For any canonical variable F, the Poisson brackets {F, H[N}} and {F,H[N]} are equal 
modulo the constraints H a (x) = = Hk(Q. Thus, the smeared constraint (3.33) is the 
canonical generator of DiffE on the constraint surface H a (x) = = Hk{(). 

D. Noether Charges 

The Noether charges Q[<fi] and Q[0] from Section 2b can be expressed in terms of the 
E-variables or the iS-variables: 

Q[4>] = ^d 3 xcf>(Z(x))U v (x) = ! 5 d 3 C0(C)iyC) , (3.34a) 

Q[6] = f d 3 x6(Z{x))n#{x) = f d 3 C0{C)U4() . (3.34b) 

JE JS 

Q[<f>] and Q[6] generate canonical transformations that yield changes in the dynamical vari- 
ables due to deformations of the hypersurfaces of constant (p and $, respectively. Among 
the canonical S- variables, only cp is affected by Q[(f>] and only i? is affected by Q[6]. 
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The Noether charge 

Q[i} = - f d 3 xe(Z(x))P k (x) (3.35) 

from Eq. (2.18) can be written in terms of the 5-variables as follows. First, write the 
momentum constraint as 

H a (x) = P k (x)Z\ a (x) + H+(x) , (3.36) 

where H£(x) contains the contributions from the S tensor fields tp, U^, •&, U$, g a b, and p ab . 
This expression can be solved for Pk(x), 

P k (X(0) = H a (X(C))X a M - H+(X(())X a 4() , (3.37) 

and the result inserted into Eq. (3.35) along with a change of integration variables x a = 
X a ((). The first term can be rewritten using Eqs. (3.15) and (3.23), while the second term 
involves a straightforward mapping of S-tensor fields to the fluid space. The result, 

Q[£] = [ d 3 (e(0(Pa(OX a ,k(0 + H+(0) , (3.38) 

can be recognized as the canonical generator of fluid space diffeomorphisms, DifLS. That is, 
for any canonical variable F, the Poisson bracket {F, Q[$]} gives the change in F due to an 
infinitesimal diffeomorphism of S generated by the vector field £(£)■ If F is one of the fluid 
space tensor fields then the Poisson bracket {F, Q[£\} equals the Lie derivative £?F. For 
the variables X a (() and P a ((), the transformation generated by Q[$] yields 

{X a ((),Q[Z}} = e(()X a , k (0, (3.39a) 

{Pa((),Q[£]} = (e(C)Pa(0), ■ (3.39b) 



Thus, X a (() transforms as a set of scalar fields on S and P a (C) transforms as a set of scalar 
densities on S. 



IV. PERFECT FLUID ACTION S 

A. Routh's Procedure 

Given an equation of state p(n, s) or p(/i, s) and the corresponding definitions (2.15), the 
Hamiltonian constraint H±(() of Eq. (3.18) depends explicitly on the canonical iS-variables 
9ke(C)j ft ke (Oi anc i n^(C), and also on s, H%, and nor/i (depending on the equation of 
state used). The variable n or \x is determined as a function of H k , s, and the canonical S- 
variables g k g(C) anc i ^l<p(C) through the auxiliary equation (3.19). The entropy per particle 
s is defined by s := II^(^)/II ¥ ,(^), and H% is determined as a function of the canonical «S- 
variables X a ((), P a (C), 9u(C)j anc i ^^(C) according to Eq. (3.20). Thus, the Hamiltonian 
constraint H±(Q depends on the canonical iS-variables g k ([C)i ^ ke (C)j X a ((), P a (C)i n^(^), 
and II^(C). It does not depend on <p(() or #(£)• Likewise, the momentum constraint 
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Hk{() of Eq. (3.23) does not depend on <p(C) or #(C). It follows that the Hamiltonian 
J s d 3 ((N ± (()H ± (() + N k (()H k (C)) for the combined perfect fluid/gravity system does 
not depend on (f(C) or i?(C); these variables are therefore ignorable (cyclic). 

Since 93(C) and ?9(C) are ignorable, they and their canonical conjugates can be eliminated 
as dynamical variables through the standard procedure due to Routh ||. Accordingly, the 
equations of motion for ip(C) an d #(C) imply that Tl v (Q and II#(C) are constant in time. 
Then 11^(0 and n^(C) can be interpreted as fixed, prescribed fluid space densities, and the 
total time derivative terms Tl v ip and Tl^fi can be discarded from the action (3.22). The 
resulting action is 

S[X a > P 8 ,0 M ,p M ,JV J -,iV*] 

= l m dtJ s d 3 c{Pa(0X a (0+P ke (09 k i(0 - N\()H ± (() - N k (()H k (0 } . (4.1) 

The Hamiltonian and momentum constraints that appear here are defined by the same 
expressions (3.18), (3.23) as before, but now H v (() and n^(C) are viewed as fixed densities 
on the fluid space S. 

The symmetries (2.12b) and (2.12c) of the S action are absent from the new action S, 
and the corresponding charges (3.34) play no role. On the other hand, the action S, like 
S, is invariant under arbitrary diffeomorphisms of the fluid space S as long as II ¥ ,(C) and 
II^(C) are transformed as densities. In this sense the description of perfect fluids provided 
by S is independent of the choice of coordinates ( k on S. However, this invariance involves 
a transformation of the fixed quantities 11^(0 and 11,? (() in addition to a transformation 
of the dynamical variables. Therefore it does not, in general, correspond to a conserved 
Noether charge. Only the subset of fluid space diffeomorphisms £ that satisfy 

£fl v = , £fl$ = (4.2) 

give rise to a conserved charge 0. This charge is given by expression (3.38) with the terms 
proportional to 11^(0 and 11^ (() dropped: 

Q[£\ = f d 3 ce(o(Pa(c)x a , k (o+m(o) ■ (4.3) 

This charge generates a symmetry of S through the Poisson brackets for any fluid space 
vector £ that leaves the particle number density 11^(0 and entropy density n#(£) invariant 
under Lie transport, Eq. (4.2). 

B. S Variables for S 

Now apply the point canonical transformation 

(X"(C), g M (0) ^ (Z k (x), g ab (x)) (4.4) 



to the action (4.1). For the coordinates X a ((), flf^(C), this is the inverse of the point 
transformation (3.7) from the canonical S-variables to the canonical S- variables. The term 
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P a {C)X a (() in the action can be rewritten using Eq. (3.10). Application of the relationship 
(3.12) (derived in Appendix B) to the gravitational kinetic term leads to 



, d i (p ke (()g ke (0= / d A x (p ab {x)g ab {x)-H G a {x)X\ k {Z{x))Z k {x 
is Js v 

Together these results yield 

d 3 ((Pa(C)X a (0+P M (09U0 



ryk 



d i x(P k (x)Z k (x)+p ab (x)g ab (x) 



(4.5) 



(4.6) 



where the momentum conjugate to Z k [x) is 

dZ(x 



Pk(x) 



dx 



P a (Z(x))X a , k (Z(x)) + H G k (Z(x))) . 



(4.7) 



Note that the momentum Pk{x) defined here differs from the momentum (which was given 
the same notation P k {x)) of Section 2b. This is clear from the simple observation that 
Eq. (3.16), which defines the relationship between P a {() and the "old" Pk{x), contains the 
variables (p, $, and their conjugates. These variables do not appear in the relationship 
(4.7) between P a (() and the "new" Pk{x). Therefore the new canonical variables Z h (x), 
Pk{x), g a b(x), and p ab (x) are not simply a subset of the canonical S-variables that appear 
in the Hamiltonian version (3.1) of the velocity-potential action S. Nevertheless, the set 
Z k (x), Pk{x), g a b{x), and p ab (x) defined here will be referred to as S-variables in order to 
distinguish them from the iS-variables X a ((), P a (C)) 9ki(()i anc ^ P H (0- 

The canonical transformation (4.4) can be accompanied by a change of variables (3.21) 
in which the Lagrange multipliers N (£) and N (£) are replaced by N ± (x) and N a (x). 
Then the action (4.1) becomes 



S[Z fc ,P fe ,<7 o6 ,p a, \A^,7V' 



dt / <Px{ P k (x)Z k (x) +p ab (x)g ab (x) - N ± (x)H ± (x) - N a (x)H a (x) 



where the momentum constraint is 

dZ(x 



HJx) 



dx 



H k (Z(x))Z k , a (x) = P k (x)Z k , a (x) + H°(x) 



The perfect fluid contribution 

dZ(x) 



H:(x) 



dx 



H F k {Z{x))Z k A {x) = P k {x)Z h 



[x 



(U 



(4.9) 



(4.10) 



differs from the previous expression (2.14a) in that the terms involving the coordinates </?, 
$, and their conjugates are absent. 

The Hamiltonian constraint in the action (4.8) can be written in terms of the E-variables 

as 



HAx):-- 



dZ(x) 



dx 



H ± (Z(x)) = H°{x) + J^m + H^H[ - ^gp 



(4.11) 



Here, H^x) is a function of Z(x) defined by 

dZ(x) 



Hp(x) 



dx 



n^Z(x)) , (4.12) 



with 11^(0 a fixed density on S. In H±(x), \x and p are determined as functions of s and 
either n or \x by the equation of state, either p(n, s) or p(fi, s), and the associated definitions 
(2.15). The entropy per particle s is a function of Z(x), s = s(Z(x)), as determined by the 
fixed scalar function 

.(0==§| (-3) 

on the fluid space S. The auxiliary equation 

TT H F n ab R F 

^ = n.l + ^f-^ (4.14) 



determines either n or //, depending on the equation of state. The auxiliary equation here 
is formally identical to Eq. (2.16), but now H£(x) and IT^ (x) are defined by Eqs. (4.10) and 
(4.12), respectively. Thus, n or /i is determined as a function of g a b(x), Pk{x) (which appears 
in H£(x)), and Z k {x) (which appears in Il (/ ,(x), s(Z(x)), and H£(x)). 

With the help of Eq. (4.7) the conserved charge (4.3) can be written in terms of the 
E-variables as 

Q[£] = - j^xe(Z{x))P k (x) . (4.15) 

— * — * 

Again, £ must satisfy the conditions (4.2). The symmetry generated by Q[£\ with such a 
vector £ corresponds to the transformations 5Z k = — £ (Z) and SPk = Pi£, \SZ\ 

C. Lagrangian S Action 

The Lagrangian form of the action (4.8) is obtained by eliminating the momenta Pk(x) 
and p ab (x). For the gravitational variables, elimination of p ab (x) yields the usual Einstein- 
Hilbert action. For the perfect fluid variables, variation of the action with respect to Pk(x) 
gives the equation of motion 

= Z* - JV* ,-N^g LP- - g- + nH ^fA , (4.16) 

vy \ 3P k dP k A.nj ) { ' 

where the auxiliary equation (4.14) has been used to eliminate the expression H^g ab H^. 
By applying the relationships (2.15) and using the fact that s = s(Z(x)) is independent of 
Pk{x), one finds that the two terms involving dfi/dPk and dp/dPk cancel one another. Thus, 
the Pk{x) equation of motion becomes 



Z k - N a Z k A = N L ^n ay ' b . (4.17) 
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With this result, the perfect fluid contribution to the action (4.8) reduces to 

S F [Z k -g ah ,N L ,N a ] = f dtf cPxN-Ly/gip-niJ,) . (4.18) 

Here, n, p, and p are determined by the equation of state, the relationships (2.15), and the 
solution of the auxiliary equation (4.14). In the auxiliary equation, Pk{x) is replaced by the 
solution of the equation of motion (4.17). 

The action (4.18) can be put in covariant form by interpreting t and x a as spacetime 
coordinates y° = t, y a = x a . The spacetime metric ^^ is constructed from the spatial 
metric, lapse, and shift according to the usual ADM decomposition, 

l a pdy a dyP = -{N^fdt 2 + g ab {N a dt + dx a ){N b dt + dx b ) . (4.19) 

This implies N^^fg = ^—7. 

If the equation of state is specified by the function p(n, s), the relationships (2.15) show 
that the covariant action S F becomes 

S F [Z h ; la p] = -Jd 4 yV^?p(n,s) . (4.20) 

Here, the entropy per particle is a given function of the fluid Lagrangian coordinates, s = 
s(Z(y)). In order to provide for a covariant specification of n, first rewrite Eq. (4.12) as 

^ = ^e abc Z\ a Z\ h Z\jl v (Z) 

= ^yg^ bc Vkt m (Z)Z k , a Z e , b Z m , c , (4.21) 

where e abc is obtained by raising indices on the volume form on E. Also, r^ m are the 
components of a three-form 

V = ^Vke m dZ k A dZ e A dZ m = U^d 3 Z (4.22) 

on the fluid space S, so that 77123 (Z) = II ¥ ,(Z) is the fixed density on S. Now consider the 
spacetime vector density 

J a := -^^e a ^ p r 1Um Z k ^Z\ a Z m , p , (4.23) 

where e a/3crp is obtained by raising indices on the spacetime volume form. Note that J a 
is tangent to the fluid flow lines, since it is orthogonal to the gradients of the Lagrangian 
coordinates. The t component of J a is J' = n^, while the spatial components are 

J a = -^e abc Vkim Z k Z\ b Z m ,c ■ (4.24) 

Now, the equation of motion (4.17) for Pk{x) can be used to obtain the result 

r + N ^ = -(^^)^K. (4.25) 

V p J J* 
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Using the ADM split (4.19), the squared magnitude of J a becomes 

\J\ 2 ■= -J a 7af3 jP = (A^J 4 ) 2 - g ah {J a + N a J l ){J b + A 6 J*) 

= (n^) 2 (4.26) 

where, again, the auxiliary equation (4.14) has been used. Thus, J a is the densitized particle 
number flux vector, J a = ny/—^U a . In the action (4.20), n is the function of Z k defined 
by n := \ J\j\J— 7, with J a defined by Eq. (4.23). The action (4.20) is discussed in Ref. 0, 
where it is confirmed that the variation 5S F leads to the correct perfect fluid stress-energy- 
momentum tensor and equations of motion. 

The Lagrangian action S, like its canonical counterpart, is invariant under fluid space 
diffeomorphisms that leave the particle number density and entropy density invariant. That 
is, the action S is invariant under transformations SZ k = — £ (Z) that satisfy 

£^q = , £%s = . (4.27) 

Indeed, with the definitions of n and s as functions of Z k , it follows that || 

Sn = ~U a e a ^{£ t n Um )Z\pZ\ a Z m , p , (4.28a) 

6s = -£%s . (4.28b) 

Thus, Sn and 5s, and in turn SS, vanish for £ satisfying Eq. (4.27). 

If the equation of state is specified by the function p(//, s), then the integrand in the 
action (4.18) can be expressed as the product of y/— 7 and p(//, s) — /xn(/x, s), where 
n(fj,,s) := dp(fi,s)/dfi. The chemical potential \i can be defined by the solution of the 
equation n(/j>, s) = \ J\/y/— 7. This gives /z(n, s), where n is now viewed as shorthand nota- 
tion for I J\/y/— 7. By substituting /i(n, s) for fi, one finds that the integrand reduces to the 
product of y/— 7 and p(/J*(n, s),s) — /i(n, s) n(fj,(n, s),s) = —p(n, s), and the resulting action 
simply coincides with the action (4.20) described above. Note that the process of solving 
for fi(n, s) and substituting the result back into the integrand is equivalent to performing a 
Legendre transformation in which the equation of state is changed from p(fi, s) to p(n, s). 
Thus, we reach the obvious conclusion that a Lagrangian S*-type action for a perfect fluid 
with a given equation of state p(p, s) can be obtained by rewriting the equation of state as 
a function p(n, s) and then using the action (4.20). 

D. Hybrid Action S F [Z k ,fr,-y af3 } 

It is possible, of course, to eliminate just one of the ignorable coordinates, either ip or 
i9, from the action S while retaining the other as a dynamical variable. This leads to a 
hybrid action that is, loosely speaking, half way between the velocity-potential type action 
S and the action S discussed above. In this subsection I will treat the case in which <p and 
its conjugate are eliminated, and i9 and its conjugate are retained. This leads to an action 
of the S type in the sense that the fluid four-velocity is expressed in terms of the totally 
antisymmetric product of gradients of Z k , as in Eq. (1.3). 

21 



The variables V'(C) an d ^1<p(C) are eliminated from the action (3.22) by dropping the 
term n^y? and elsewhere interpreting 11^(0 as a fixed density on the fluid space S. The 
resulting action can be expressed in terms of E-variables by applying the point canonical 
transformation 

(X a (0, 0(C), g ke (()) -> (Z k (x), 0(x), g ab (x)) . (4.29) 

The momenta conjugate to ~&(x), g a b(x), an d Z k (x) are IL#(x), p ab (x), and 

dZ(x) 



Pk(x) 



dx 



{P a {Z{x))X\ k {Z{x)) + U4{Z{x))Q, k {Z{x)) + H G k (Z(x))) , (4.30) 



respectively. Associated with the Lagrange multiplier N a (x) is the momentum constraint 
H a (x) = H^{x) + H°{x) with the perfect fluid contribution 

H*{x) = P k {x)Z\ a {x) + n*(a;)# a (s) . (4.31) 

The Hamiltonian constraint associated with the Lagrange multiplier N ± (x) is given by 
Eqs. (4.11), (4.12), and (4.14). In this case the entropy per particle, 

s = TW ' ( } 

is a function of IL#(x) and Z k {x) (through U v (x)). 

The symmetries of this action are generated through the Poisson brackets by the con- 
served charges 

Q[0] = ( d 3 xTl4x)6(Z(x)) , (4.33a) 

JY, 

Q[l] = Jj 3 xP k (x)e(Z(x)) , (4.33b) 

— * 

where £ leaves the fluid space number density invariant, £^11^ = 0. The charge (4.33a) 
generates a deformation of the constant d surfaces, 8d = 0(Z), while the charge (4.33b) 
generates a fluid space diffeomorphism, 5Z k = —£ k (Z). 

The action can be written in Lagrangian form by eliminating the momenta P k (x), H#(x), 
and p ab (x). Elimination of p ab (x) leads to the usual Einstein-Hilbert action for the grav- 
itational field. Varying the action with respect to P k (x) produces the equation of motion 
(4.16), which reduces to Eq. (4.17). With this result, the perfect fluid contribution to the 
action becomes 

S F [Z k 1 d^ d -g abl N ± 1 N a } 

= j^ dt jf d 8 *! tfVff (P - W) + n, (^ - N a d A - N ± V9 nH "j£*' b ) } ■ (4-34) 

At this stage it is convenient to make a change of variables in which s = ILi/Tl^ replaces 
Iltf as an independent variable. (Thus, for the moment, the momentum ILj will not be 
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eliminated from the action, but will be hidden in the new variable s.) Assuming now an 
equation of state of the form p(n, s), the action is 

S*[Z k ,&,s;g ab ,N ± ,N a ] 

= J^-N^pin, s) + sU v (d - N^ >a - N ± V9 nH °j£*' h ) } ■ (4-35) 

Again, define the vector density J a of Eq. (4.23). With the result (4.25) it is straightforward 
to verify that 



_ n HFgOt>#. 



J a # a = Il ip [#-N a #, a -N ± ^g * y '" ) , (4.36) 



so that the action simplifies to 

S F [Z k ,0,s;j O 0]=Jd*y( r V z: 7P(n,s)+sJ a tOt ) . (4.37) 

Here, n = \ J\/ \/— 7 depends on Z k through the definition (4.23) for J a . 

The entropy per particle s (and hence the momentum ILj) can be eliminated from the 
action (4.37) by use of the solution of the s equation of motion, namely, 

= -^^ + ™,«- (4.38) 

From the first law Eq. (2.1), we have dp/ds = nT so that d is related to temperature by 
T = $ ta U a (Eq. (2.5)). It follows that the solution of the s equation of motion (4.38) has the 
form s(n,T), the entropy per particle as a function of number density n and temperature 
T. If this result is substituted into the action (4.37), the integrand becomes the product of 
^—7 and a function of n and T, 

n a(n, T) := nT s(n, T) — p(n, s(n, T)) . (4.39) 

The function a(n,T) is the physical free energy of the fluid. The action (4.37) therefore 
reduces to 

iF[ ryk 



S F [Z\& lla p\ = - / d*y^na(n } T) , (4.40) 



where n = \J\/y/— 7 and T = ^, Q J Q /| J\ with J a defined as a function of Z k by Eq. (4.23). 
This action is invariant under the symmetry transformations 5$ = 9(Z) and 5Z k = — £ (Z), 
where £ leaves the fluid space three-form r\ invariant under Lie transport, £p] = . 

The action (4.40) is the one discussed by Kijowski et al. [0. The equation of state 
is specified by giving the physical free energy a as a function of number density n and 
temperature T. 
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APPENDIX A: HAMILTONIAN ACTION WITH EQUATION OF STATE P(p, s) 

Let t denote a scalar function on spacetime whose gradient is nonzero and (future point- 
ing) timelike. The t = constant surfaces foliate spacetime into spacelike hypersurfaces with 
unit normal n a = —N^t^, where N 1 - := (—t )a 'j al3 t^)~ 1 ^ 2 defines the lapse function. Also in- 
troduce the time vector field t a such that t^ a t a = 1 and the projection tensor g^ := 5p+n a np 
onto the leaves of the foliation. With the shift vector defined by N a := g^t 13 , the time vector 
field can be written as t a = N- L n a + N a . 

Using the above relationships, the vector V a can be decomposed as 

V a = (g% - n a np)VP = (gV) a + (t a - N a )V l , (A.l) 

where the definitions (gV) a := g'SV^ and V 1 := t j0 ,V a are used. It follows that the chemical 
potential \V\ can be written as 

\V\ := spV^V? = ^(iW) 2 - (gV) a (gV)<* . (A.2) 

The action (2.9) then becomes 

d 4 yV=l{-p(\V\,s) + U ^ S \ vH a + (gVr - V l N a ){p^ + s#, a - W k Z\ a )} , 

(A.3) 



S 



where the definitions n(p, s) := dp(p, s)/dp and p(p, s) := pn(p, s) —p(p, s) have been used. 
Now tie the spacetime coordinates to the foliation by choosing t as the time coordinate 
and choosing spatial coordinates such that d/dt is the time vector field t a . The t-components 
of (gV) a and N a vanish, and the action becomes 



S F [{gV) a X,ip,$,s,W k ,Z k -N L ,N\ 

= J^dt J^xN-Ly/g^-pdVls) + W(l ,^j S V (y + stf - W k Z k ) 

+ W(l ,^| S) ((^) a - V l N a )(p A + sd, a - W k Z\ a )} . (A.4) 

The form of this expression suggests a change of variables in which V 1 , s, and W k are replaced 

by 

IT, := N^^^-V* , (A.5a) 

II* := N^^^V's , (A.5b) 

P k ■= -N^^^-V'W, . (A.5c) 

These variables will become the momenta conjugate to tp, i9, and Z k , respectively. 
Now introduce the variables 
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C ■= {j ^r- , (A.6) 

which can be viewed as replacements for (gV) a . In terms of the new variables (A. 5), (A.6), 
the action becomes 

S F [e, V, K v , #, n*, Z k , P k ; N\ N a , g ah ] 

= J^dtJj\{ll^ + Uj + P k Z k -N^^ P {\VlTl»/n^ + {e-N a )H F a } , (A.7) 

where s = U#/U v (as follows from Eqs. (A. 5a) and (A. 5b)) and H F is the fluid contribution 
to the momentum constraint, 

H F := n^,a + n*#, a + PkZ k , a . (A.8) 

In the action (A.7) \V\ is determined as follows. Definition (A.6) can be combined with 
Eq. (A.2) to yield 

V* = , ' ' . (A.9) 

v / (A r± )2 _ ^" 

Then the momentum (A. 5a) becomes 

, «VM. (A.10) 



This can be viewed as an equation that determines \V\ as a function of IT^, H#, £ a , N- 1 , and 

<7a6- 

The next step is to eliminate £ a from the action (A.7) by substituting the solution of the 
£ a equations of motion. The £ a equations of motion are found by varying S with respect 
to £ a , keeping the other variables fixed. For this calculation, the variation of |V| is needed. 
From Eq. (A. 10) this is determined to be 

fti(|y|,iwn y ) MiniWn y ) \ , A . 

Using this result and the relationship dp(fi,s)/dfi = fidn(fi,s)/dfi in the variation of the 
action (A.7), the £ a equations of motion are obtained: 

I VI H F 

1 ' & = -jf- . (A.12) 



Here, \V\ is considered to be a function of IT^, n#, £ a , N- 1 -, and g a b through the solution of 
Eq. (A. 10). Equation (A.12) can be squared and, after a bit of algebra, written as 



N A 
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H F g ab H F . A . 

ay " (A.13) 



Then Eq. (A. 10) becomes 



^ = //(|r|.ii,/n 



H '*** H i (A.14) 



V~9 ""'""*'* *»\ \v\mi ■ 

This equation determines \V\ as a function of H^, 11$, H£, and g a \,. 

The implicit solution (A. 12) of the £ a equations of motion can be inserted into the action 
(A. 7) to eliminate the variable C, a . Note that C, a only appears in the fluid contribution to 
the Hamiltonian constraint; this is 

Hi = Jg\V\n(\V\,ILe/Ily) - ^f- - ^p(\V\,U$/U v ) , (A.15) 

where p(fi, s) = fin(fi, s) — p(/i, s) has been used. The term involving the number density 
n can be rewritten using Eq. (A. 10), and the term involving C, a H^ can be rewritten by 
contracting Eq. (A. 12) with £ a . This leads to the expression 

H * = n ^l 77 Z l - V9P(\V\,Tl,/Tl^ , (A.16) 

y( N ) ~~S ?6 

which, when combined with Eq. (A. 13), becomes 



H* = J\V\*IP v + HZg<*H?- y /gp{\V\,Ile/Il <f> ) . (A.17) 

Here, \V\ is determined as a function of 11^, 11^, H£, and g a b through Eq. (A.14). 

With the result (A.17) the action (A. 7) is brought to standard Hamiltonian form with 
contributions (A.17) and (A. 8) to the Hamiltonian and momentum constraints. The chem- 
ical potential \V\ that appears in Hf_ is determined by the auxiliary equation (A.14). This 
agrees with the action described in Eqs. (2.13)-(2.17) (for the equation of state p(n,s)) of 
the main text. 



APPENDIX B: DERIVATION OF EQUATION (3.12) 

In this appendix the parameter dependence of tensors and mappings will be denoted 
by subscripts. Thus, for example, the one-parameter family of mappings from S to E is 
denoted X t , and the generic one-parameter family of tensors on E is denoted Xt- 

Let V t denote the vector field on E defined in terms of local coordinates by 



K( X ) := (|^(C)) 



(B.l) 

C=Zt(x) 



For each value of t, the integral curves of V t define a one-parameter a family of diffeomor- 
phisms Mt,o : E ^ E. By definition J\ft,a satisfies the differential equation 

^ = VM„) , (B.2) 
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where Mt,o = It, is the identity mapping on E. By combining the previous two equations, 
we obtain 



d 

da 



Krl 



X 



d 



xm 



C,=Zt{x) 



d 

da 



(7=0 v dt 

It readily follows that up to order a 2 , J\ft,a is given by 

Mt, a = X t+a o Z t + 0(a 2 



X^(Zt(x)) 



a=0 



(B.3) 



(B.4) 



The 0(a 2 ) corrections to Ht,a are not needed in the analysis that follows. 

Consider first the case in which Xt is a one-parameter family of scalar fields on E. The 
Lie derivative of \t with respect to V t is given by 



( £ v t Xt) 



d 

da 



Xt(M,a(x)) 



(7=0 



d 

da 



Xt(Xt + a(C)) 



a=0, C=Z t (x) 



(B.5) 



In order to simplify this expression, denote the pullback of Xt{ x ) to the fluid space S by 
Xt(() '■= Xt(X t (()) and compute 

di Xt{C) = da' Xt+ ^ C) 



cr=0 



d 



da 

d_ 

dl 



Xt + a(X t+(T (()) 



a=0 



Xt{x) 



d 



+ ^Xt(X t+a (Q) 
x=x t (0 da 



a=0 



Thus, the Lie derivative (B.5) equals 



( £ y^)\ r -j t xt(x) + (j t x t (0 



C=Zt(x) 



(B.6) 



(B.7) 



This result expresses the Lie derivative of Xt along Vt in terms of the derivatives of Xt and 
Xt with respect to the parameter t. 

The result (B.7) can be extended to arbitrary one-parameter families of tensor fields Xt- 
For example, consider Xt(%) to be a covariant vector and compute 



(£v t x> 



d_ 

da 



d (x>. 



a=0 



= z* t — [x* +ff xt{x t+a {0) 

Here, Eq. (B.4) has been used to express 

Nl a = z; o x; +a + o(a 2 

A calculation similar to that in Eq. (B.6) yields 



ff =o, (=Z t (x) 



d 
It 



*(o =*?(!*(*)) 



x=X t (0 



+ ^(x; +ffX t(x t+(T (0) 



a=0 



(B.8) 



(B.9) 



^B.10) 
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for the fluid space covector Xt(0 := XtXt(Xt(C))- Therefore the Lie derivative of Xt is 



£v t Xi 



-jM x) + <jMQ 



[B.ll] 



C,=Zt{x) 

Likewise, the Lie derivative of a contravariant vector is given by 

.// d 



£v t x-i 



d_ 

~dt 



x,W + x;(- x «(o) 



C=Zt{x) 



(B.12) 



where X' t is the derivative mapping of contravariant vectors. 

Armed with the results (B.7), (B.ll), (B.12) and their generalizations to arbitrary ten- 
sors, Eq. (3.12) is obtained as follows. For defmiteness, consider the case in which Xt{%) 
is a one-parameter family of covariant vectors and ir t (x) denotes the canonically conjugate 
one-parameter family of contravariant vector densities on E. Contract Eq. (B.ll) with 7it{x) 
and integrate over E to obtain 



d 3 xn t (x)—xt(x) = j^d 3 x(n t (x)Z;(—Xt(C) 



C=Zt(x) 



C=Zt(x) 



TT t (x)(£ Vt Xt 



V t a (x)H a (x;7r t ,Xt]) ■ (B.13) 



Here, H a (x;n t ,Xt\ is the contribution from Xt an d ftt to the momentum constraint. (Inte- 
gration by parts is used to obtain the second line of Eq. (B.13) from the first. It is assumed 
that the space manifold E is compact without boundary, so that no boundary terms ap- 
pear.) Now rewrite the right-hand side of Eq. (B.13) as an integral over the fluid space S 
by changing integration variables, x a = X®((). This gives 

J^xn{x)j t xt{x) = J s d 3 ((ir t (Oj t Xt(() - V k t (OH k ((;ir t ,X t }) ■ (B.14) 



The vector field Vt{() is defined in terms of local coordinates by 



V*(():= (V t a (x)Zl a (x, 



d 

x=X t (0 \di 



** a (C)) (%*(*)) 



(B.15) 



x=X t (C) 



where the definition (B.l) has been used. Inserting this expression for Vt{C) ^° Bq. (B.14) 
yields Eq. (3.12) from the main text. The derivation for tensors of arbitrary rank follows 
along similar lines. 
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